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Analytic results are presented for preheating in both flat and open models of chaotic inflation, 
for the case of massless inflaton decay into further inflaton quanta. It is demonstrated that preheat- 
ing in both these cases closely resembles that in Minkowski spacetime. Furthermore, quantitative 
differences between preheating in spatially-flat and open models of inflation remain of order 10~ 2 
for the chaotic inflation initial conditions considered here.n 
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I. INTRODUCTION 



Recently, a new view of the post-inflationary reheating period has been established. [|l|j2| In place of the original 

view, in which the inflaton decayed perturbatively J3[, an inherently nonperturbative, highly-efficient resonance has 

been investigated. The new theory of reheating now involves three distinct stages: an oscillating inflaton sets up 

a parametric resonance in its decay to some boson species; this explosive stage has been termed "preheating." Q| 

— i ', (See also |^,[I| Q.) Next, these far-from-equilibrium decay products themselves interact and decay, as can be stud- 

►> ied along the methods in |H|. And finally, the decay products thermalize, completing the reheating process. Such 

■^j- \ explosive preheating can radically change the thermal history of the universe following inflation; some non-standard 

effects associated with the new preheating picture include the possibility for non-thermal symmetry restoration [15| , 

supersymmetry breaking and GUT-scale baryogenesis fll6| , and the amplification of gravitational radiation following 

the preheating period Q . 

Because preheating is a non-linear, non-equilibrium process, particle production in most specific models must be 
analyzed numerically, and several results have been reported for both Minkowski and spatially-flat expanding space- 
times. P|,n,P, |l0|Jlq 1 Much of the analytical literature to date has treated the expanding background case by means of 
i-Q two approximations: (1) that the oscillating zero-mode of the inflaton oscillates as a cosine function, ip ex cos(mi), 

*tH with m the mass of the inflaton, and (2) that the expansion may be neglected for frequencies m 3> H, where H is 

O i' the Hubble parameter. When these approximations are made, the equation of motion for the quantum fluctuations 
(associated with particle production) can be cast in the form of a Mathieu equation. Solutions of the Mathieu equation 
generically exhibit an infinite hierarchy of resonance bands; for wavenumbers k within these resonance bands, the 
solutions grow exponentially in time, driving the explosive, resonant particle production. As discussed in J9|, however, 
in the context of Minkowski spacetime, the first of these approximations can lead to large quantitative errors when a 
quartic self-coupling exists for the inflaton. In these cases, the zero-mode evolves as an elliptic function in time, and 

jrt ■ solutions for the fluctuating fields reveal only one single resonance band in k. 

In this paper, we extend this analytical study to the case of inflaton decay into inflaton bosons in an expand- 
ing Friedmann-Robertson- Walker spacetime. For a quartic interaction potential and a massless inflaton, the time- 
evolution of the entire system (zero-mode, quantum fluctuations, and background spacetime) can be studied con- 
sistently with analytical methods. Because in this case the observed spectrum of cosmic microwave background 
anisotropies places strict limits on the self-coupling, A ~ 10~ 12 , we may study this nonperturbative, non-equilibrium 
system by means of analytic approximations to the full, non-linear equations of motion. In ||, the large- A^ approxi- 
mation is employed to study preheating analytically in Minkowski spacetime. Here, we rely on a Hartree factorization 
to study preheating in an expanding spacetime. As demonstrated explicitly in [|19|| , these two approximation schemes 
are closely related, and in the case of preheating amount only to the substitution A — ► 3A in the equation of motion for 
the fluctuating field. (See also p0[ .) Numerical results have indicated that the Hartree approximation can break down 
for quantitative results far from the weak-coupling range (that is, for inflaton decay into a distinct species of boson, 
with coupling constant g> A) m. However, for the case studied here, of inflaton decay into other inflaton quanta, 
the Hartree approximation should remain quantitatively reliable. With this analytical framework, it is also easy to 
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understand numerical results ]7|,p8| which indicate that no parametric resonance can occur for a massive inflaton in 
expanding spacetime, if its only decay channel is into inflaton quanta. 

We present analytic results for preheating for both ordinary chaotic inflation and for chaotic inflation within a 
model of open inflation. |El| The growth of the backreaction from the created quanta on the oscillating zero-mode is 
calculated, as well as the maximum number of quanta produced during preheating. As discussed below, for preheat- 
ing of a massless inflaton into massless inflaton quanta, the Ricci curvature scalar vanishes during preheating. Thus, 
the preheating dynamics for this model in an expanding spacetime remain conformally equivalent to the Minkowski 
spacetime case; it is not surprising, therefore, that the results presented here share the same form as many of the 
Minkowski results studied in H. The change from large- N to Hartree methods, however, does change some of the 
details of the quantitative analysis. Also, as demonstrated below, the spectrum of produced quanta in the open 
inflation scenario takes the same form as that for quanta produced in an expanding, spatially-flat model; for chaotic 
inflation initial conditions, there are only very small numerical deviations between the two situations. 

In section II, we present the specific model to be analyzed here, discuss the approximation scheme to track the 
quantum backreaction, and choose an appropriate initial vacuum state with reference to which we can measure the 
resonant particle production. Section III includes analysis of the background spacetime dynamics during the pre- 
heating phase. In section IV, we solve consistently for both the oscillating zero-mode of the inflaton field and the 
coupled quantum fluctuations, for both the flat and open universe cases. In section V, we compare numerically the 
preheating spectra for flat and open inflation, and provide further consistency checks on some of the approximations 
made throughout the analysis. Concluding remarks follow in section VI. 



II. DYNAMICS OF THE MODEL 

We only consider the case of inflaton decay into inflaton bosons, due to a non-linear self-coupling. The Lagrangian 
density thus may be written: 




L iR-lg^d^d^-Vi® 
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We will consider the case of an additional non-minimal coupling between the inflaton and R(t), the Ricci curvature 
scalar, in section III. The line element for a general Friedmann-Robertson- Walker (FRW) spacetime may be written 

ds 2 = -dt 2 + a 2 {t)h io dx l dx\ (2) 

with 

h ij dx i dx j = dx 2 (K = 0), 

= dr 2 + sinh 2 r (d9 2 + sin 2 Odcf) 2 ) (K = -I), (3) 

and K = for a flat universe, K = — 1 for an open universe. (Naturally, the angular coordinate <f> should not be 
confused with the inflaton field 4>(t, x).) In this metric, the equations of motion follow: 
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where H = a/a is the Hubble parameter, V 2 is the comoving Laplacian operator, and dots represent derivatives with 
respect to cosmic time, t. 

Making the familiar decomposition between the "classical" and fluctuating portions of the inflaton field, 

0(t,x)= ¥ »(t) + *0(t,x) J (5) 

we see that the potential of equation (||) contains a coupling between tp and 8<p. It is this coupling which, under 
certain conditions, will allow for a highly efficient transfer of energy from the "classical" to the fluctuating portions 



of the inflaton field; this transfer of energy is manifested as a rapid production of out-of-equilibrium infiaton quanta. 
Because the self-coupling strength A is constrained to be very weak for this model (A ~ 10~ 12 ), based on the observed 
anisotropies in the cosmic microwave background radiation, we may employ the Hartree factorization to study the 
non-equilibrium dynamics of these coupled systems. [|22f| This factorization amounts to making a particular choice of 
vacuum state (as discussed below), and making the following substitutions: 

W) 4 ^6(W) 2 )W) 2 ~3<(^ 2 )) 2 , (6) 

where quantities in brackets indicate expectation values of the associated quantum operators (to be defined explicitly 
below); the tadpole condition further requires {(6<fr)) = 0. 

In an open universe, there exists a physical curvature length-scale a(rf)/\K\, with K = — 1. The corresponding 
comoving curvature scale is thus simply +1. As usual, scalar fields may be expanded in eigenfunctions of the comoving 
Laplacian. Writing V 2 /p(x) = — P 2 /p(x), then in an open universe, modes with eigenvalue P 2 > 1 vary on scales 
shorter than the comoving curvature scale, and hence may be labelled "subcurvature" modes, whereas modes with 
< P 2 < 1 correspond to "supercurvature" modes. |2^,|24J As first noted in J25I], and as demonstrated below in section 
V, preheating in an open universe will only populate subcurvature modes for potentials of the form in equation (fil). 
For this reason, we may expand the fluctuating field 8(f) as follows: 



5<f>(t, x) = / dp,(k) 5<t>kjm{t) Z k jm(x), (7) 

where the eigenfunctions of the Laplacian obey 

V 2 Z kjm (x) = -(k 2 - K)Z kjm (ps). (8) 

Here k is the comoving subcurvature wavenumber, with < k 2 < oo; it is related to the ei gen value P 2 above by 
k 2 = P 2 + K, and only applies for P 2 > I. The measure for the wavenumber integral is 23 24Jj2q| 



dfl(k) = dkk 1 (K = 0), 

Jo 

= *EE (K = -l). (9) 

Ja j=o m=-j 

With these definitions, the eigenfunctions of the Laplacian may be written: 

Z kjm (x) - (2vrr 3 /V k - x (A- = 0), 

Z kjm (r, 9, <t>) = J-f[ {s 2 + k 2 y 1/2 Y jm {6, 0) 

/ -1 d\ j+1 

x (smhr) 1 -— cos(fcr) (K = -I), (10) 

\ snili r dr J 

with Yj m (Q, </>) the usual spherical harmonics. p4j2q,M 

To study particle production during the preheating epoch, we may now promote the field 5<f> to a Heisenberg 
operator: 

6<f>(t, x) = / dfi(k) 84> k (t) d kjm Z kjm (-x) + S(f>* k (t) d ] kjm Z% jm {-x) , (11) 

with d k j m and aL m the canonical time-independent annihilation and creation operators, respectively. These are 

- 2 

defined with respect to the initial Fock vacuum. With this, the expectation value (Scj) ) becomes, for both K = and 
K = -1, 
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where translational invariance allows (5(f) ) to depend only on time; we will therefore write this quantity as (8(f) (£)). 
Defining the constant 



e 2 = <<^(t )), 



(13) 



where t is the beginning of the preheating epoch, and switching to conformal time, dr) = a 1 dt, we may introduce 
the following dimensionless variables: 

r = Cr] , £ = k/C , M = m/C , 

1 Xt{r) 



VA a(T) yfC a(r) 

S(r) = [(X 2 (r)) - <X 2 (r Q ))] . 



(14) 



Then ((50 (£)) = C 2 a 2 (T)(x 2 ( r ))! an d using the factorization of equation (pi), the equations of motion for the ip an d 
X fields may be written: 



d 2 a" 

—r + a 2 (r)M 2 + iP 2 {t) + 3A« 2 (r ) + 3AE(r) 

jit a 

-^—+£ 2 -JC + cl 2 {t)M 2 - — + 3yj 2 ( T ) + 3Aa 2 (r ) + 3AS(r) 
dT z a 



V(r) = 0, 
Xi(t) = 0, 



(15) 



where primes denote derivatives with respect to r, and K, = C 2 K. It is convenient to define the frequencies Wi(t) 
and uji(t) as 



W?(t) 



AC- 



a ^ T )M 2 -—, 
lu 2 (t) = W £ 2 (r) + 3V> 2 (t) + 3Aa 2 (r ) + 3AS(r). 



(16) 



The quantity S(r) measures the backreaction of created quanta on the evolution of the oscillating zero-mode, ip( T )-> 
and S(r ) = 0. We will study solutions of the coupled equations of equation (15) in section IV, for early times when 
A£(t) may be neglected. We will also determine self-consistently in that section the time at which the approximation 
A£(t) — > breaks down. 

It remains in this section to derive an expression for the particle number operator appropriate for the non-equilibrium 
dynamics. In terms of the rescaled field Xj the quantum fluctuation operator 5<j) may be written: 



VC J a[T) 



(17) 



Xe(r) d l]m Z ijm (x) + xH T ) Ajm Z £j m ( x ) 

and, from the Lagrangian density in equation (pi), the Heisenberg operator for the conjugate field may be written: 

VCU(t, x) = f dfi(e)a(T)y/h(xj[(xi(.T) - H(t)xi{t)) a ijm Z £jm {x) + h.c] , (18) 

with T~L(t) = a' /a, and "h.c." denoting hermitian conjugate. Canonical quantization then gives the normalization 
condition for the mode functions Xi{ T ) : 



XI Xi -XeXe 



(19) 



This normalization comes from quantizing the fields 5<j) and II on a Cauchy surface. Yet for models of open inflation, 
the interior of the nucleated bubble is not a Cauchy surface for the entire de Sitter space. It has been demonstrated, 
however, that for scalar fields expanded only in subcurvature modes, quantizing as if the interior of the nucleated 
bubble were a proper Cauchy surface reproduces exactly the same result as when the fields are quantized on a proper 
Cauchy surface and analytically continued inside the bubble. [|28| For this reason, we may use the normalization in 
equation (M) for both the flat and open inflation cases. In the context of preheating, this makes sense physically, 
since we are quantizing these fields after the volume of the interior of the bubble has expanded by a factor of ~ (e 65 ) 
(according to any "observers" on the interior of the bubble), so that all causal properties of these fields should be 
specifiable with reference to the interior of the bubble alone. 

To calculate the particle number per mode, it is convenient to write 



(20) 



6<j>(t, x) = / dp,(£)5<f>i jm (t)Zi jm (x), 
n(i, x) = / dp,(£)ti.i jm (t, x)Z ijm (x). 



For K = 0, we may thus write 



1 



£jm 



(*) = 



a(r) 
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VCIl ijm (t, x) = a(r)v^(xy [(xj - H X i) a r + (xf - ft*?) a f _ £ J . (21) 

As pointed out in p4j for the K = — 1 case, we may use the purely- real representation of the Yj m {9, </>)'s, which makes 
the eigenfunctions Zij m (x.) purely real for subcurvature modes. Thus, for K = —1 we also may write expressions 
as in equation (pl|), but with the arguments of the creation and annihilation operators changed to: a^* — > a^ m and 



-i 



at 



Care must be taken when studying preheating in an expanding universe to choose an appropriate vacuum state 
with respect to which we can measure particle production. There are two independent concerns: first, preheating 
is a non-equilibrium process among interacting fields, so defining "particle" states is ambiguous even in Minkowski 



spacetimc 
spacetimc 



in addition, we have the usual ambiguity in defining "free particle" states in any non-Minkowskian 
26(| ■ At r — ► — ex) (where we take the bubble nucleation to have occurred at r = — oo), the spacetime inside 
the nucleated bubble is a de Sitter spacetime, so in the absence of interactions, the vacuum state for the 5<f> field 
should be the Bunch-Davies vacuum. If the transition from pure de Sitter expansion to the different rate of expansion 
at the time of preheating occurs adiabatically, then at the onset of preheating the \i modes should obey the initial 
conditions: xe( T o) = (2W^(t )) -1 / 2 and X^(r ) = —i(Wf(r )/2) 1 / 2 in the absence of interactions. These modes would 
represent the free-particle "adiabatic" states. p6J Furthermore, if we consider the interaction strength to be turned 
on adiabatically beginning some time before r , then these initial conditions should be replaced by: 



Xi{to) = 



1 



y/2u}i(i 



dr 



ui{t ) 



(22) 



Equation (E2) gives the initial conditions for the "adiabatic" particle states for the non-equilibrium dynamics at the 
onset of preheating. We may therefore define adiabatic creation and annihilation operators &^j m (r) and 6i\- m (T) 
through the relations: 



1 



1 



1 
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aej m {T)e J +aL m (T)e J 



VCU ijm (T, x) 
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H{t) 



/ \ —i I Ui(r)dr , , 



(23) 



The aij m (T) annihilates the time-dependent adiabatic vacuum state: d£j- m (r)|0(r) >= for all £, j, and m. These 
adiabatic creation and annihilation operators can be related to the time-independent operators &tj 



and a\ jm by 



means of a Bogolyubov transformation. As above, we replace the arguments of these operators as &ij m — > &/ and 
d]- — > a -in equation (23) for the K — case. Note that these expansions are only valid when the frequency u>i{t) 
is purely real; as demonstrated in section III, this will always be the case for the scenarios considered here. 

From the expansions in equations ( pl| ) and (23), we may solve for the adiabatic particle number per mode. The 
result is the same for both the K = and K = — 1 cases: 



^fV) = (4jm(r)&i jm (T)) 



0J£{t) 



\Xi(r)Y 



1 



u}(t) 
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(24) 



This yields the number of "adiabatic" -state quanta produced relative to the initial Fock vacuum, |0 >= |0(r o ) >. 
With this choice of vacuum state and initial conditions, the particle number for inflaton quanta does indeed vanish 
at the onset of preheating. As discussed in section III, for the models of interest here, the background spacetime will 
evolve as if it were radiation-dominated for the entire period of preheating, so that there will be no further Bogolyubov 
transformations needed to relate the particle number at the beginning to that at the end of preheating. 

With the choice of initial conditions for the \ field, equation (|2^), the quantity C 2 = (5<j) (r )) is formally quadrati- 
cally divergent, and would have to be renormalizcd with some regularization scheme. (This is essentially the zero-point 
energy divergence.) For here, we will simply note that to remain consistent, we require the energy density of the "clas- 
sical" portion of the inflaton field to exceed that of the quantum fluctuation portion at the beginning of preheating; 

- 2 

that is, </? 2 (t ) 3> (S4> {t )). This is equivalent to the requirement 



tp 2 (T ) > Aa 2 (r ), 



(25) 



using the definition of ip(r) in equation (14). We will make use of this relation below 



III. EVOLUTION OF THE BACKGROUND SPACETIME 

In this section we study the background spacetime dynamics during the period of preheating. Up until the time 
T , the fluctuations of the field Sep are in their vacuum state, and we need only consider the energy density from the 
"classical" part of the inflaton field, ip. In terms of cosmic time t, this may be written: 



1 



1 



1 



In terms of dimensionless conformal time r and the rescaled field ip, the energy density is thus 



M r ) 



c- J 



2Aa 4 (r) 



^V 4 + V'' 2 + a 2 M 2 4> 2 + Hip {Hip - 2ip') 



(26) 



(27) 



At the onset of preheating, the field ip begins to oscillate (as treated explicitly in section IV), so that averaging 
over a period of its oscillations, both ip 2 and ip will remain nearly constant. For the chaotic inflation scenario we 
consider here, then, when the mass m vanishes, the energy density at the onset of preheating will be dominated by 
the a~ A terms. Note that by working in terms of conformal time, we do not have to make any assumptions about 
the magnitude or rate of change of the Hubble parameter H(t), as we would if we worked in terms of cosmic time, t. 
Instead, we only need to use the chaotic inflation initial conditions, ip 4 (r ) 3> ip 2 (T ). We will confirm in section V 
that when m — 0, these initial conditions ensure that H <C ip during preheating. 

Thus, when m = (and only then), we may approximate the time-evolution of the background spacetime as that of 
a radiation-dominated FRW metric at the very onset of preheating. (This point is also noted in [14}.) Furthermore, 
for the case of inflaton decay into inflaton bosons, the produced quanta will also be massless when m = 0, so that over 



the entire preheating period we may keep the approximation p(r) 



(t) . Note that this behavior of the energy 



density holds even though the quanta are far from thermal equilibrium when produced. 

Using the Friedmann equation (see equation (||)), and writing p(rj) = p {a{r}) / a )~ Ag , we may solve for the behavior 
of the scale factor in terms of conformal time 77 for both the K = and K = — 1 cases. For a flat universe, the scale 
factor evolves as: 



and for K 



\l/(2fl-l) 1 
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(28) 



(29) 



In this case, by defining a — a(rj ), the Friedmann equation further implies the relation H 2 — 2/a 2 . In terms of the 
dimensionless conformal time r, setting g = 1 (for m = 0) yields the very simple result for K = 0, 1: 



a 

a 



-K. 



(30) 



This result, valid only for m = 0, means that the addition of a non-minimal coupling between the inflaton and the 
Ricci curvature scalar of the form i£i?</> 2 would not affect any of the preheating dynamics. To see this, consider the 
Ricci curvature scalar in terms of cosmic time, t: 



m = J (r) [&a + a 2 + K] 

Rewriting this in terms of the dimensionless conformal time r, 

6C 2 



R(r) 



a 2 (r) 



+ K 



(31) 



(32) 



we see from equation (|30|) that a consistent solution of the modified Friedmann equation would include R = for 
both K — and K = — 1 when m = 0. Thus, for the massless case, a non-minimal coupling would not affect the 
preheating dynamics. 

Using equation (j30|), we may further confirm that the frequency loi{t) will always remain real (see equation (P-6|)), 
which allows us to employ the adiabatic creation and annihilation operators of equation (E3). In fact, when m = 0, 
this frequency is the same for both the K — and K = — 1 cases: 



w 2 (r) = f + 3iJj 2 (t) + 3Aa 2 (r ) + 3AE(r). 



(33) 



The condition of equation ( |25| ) further allows us to neglect the 3Aa 2 (T ) term relative to the 3tjj 2 term during pre- 
heating. With these expressions for the evolution of the background spacetime, we may now study the dynamics of 
the tp and \ fields during preheating. 



IV. EVOLUTION OF THE FIELDS DURING PREHEATING 

Using equations d2q) and (BCf), we may rewrite the equation of motion for the ip field in equation fllq) as 



- t +^ 2 (t)+/C + 3AE(t) 



iJ){t) = 0. 



(34) 



In this section, we will shift the time of the onset of preheating to r = 0. With the initial conditions, ip(0) = ip and 
(dip/d,T)\ T —o — 0, the equation of motion for ip{ T ) m ay be solved in terms of a Jacobian elliptic cosine function for 
early times, when the backreaction 3AE(r) is negligible. In general, the Jacobian cosine function cn(it, k) obeys the 
differential equation E9]: 



^ + (i - 2fc 2 + ^ 2 c 2 ) 



c = 0, 



(35) 



where we have used the abbreviation c = cn(u, k). For early times, then, the ip field oscillates as 



ip(r) = i/> en ^/%l)l +SCt, 



for ip 2 > — 2/C, appropriate for a chaotic inflation scenario. 
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In terms of the time-like variable u = \/ip% -\-K.t, the equation of motion for the fluctuations becomes: 

d 2 



7 2 ^ + ^ + 3^( M) 



Xt(u) = 0, 



(37) 



for early times, when we may neglect AE(t). Here we have defined the constant 7 2 = ip 2 + /C. We may rewrite the 
Jacobian cosine function in terms of the doubly-periodic Weierstrass function p(z) by noting both that cn 2 (u, k) + 
sn 2 (w, k) = 1, and that 



sn 2 (u, k) = —^ 

k (ei - e 3 ) 



P 



u + iK'(k) 



e-z 



(38) 



where K'(k) is the complementary complete elliptic integral of the first kind, and the e, are the three constants 
associated with the Weierstrass function p(z). [E9[ These constants sum to zero: e\ + e2 + e$ — 0; in terms of them, 

2 

the modulus k may be written: k = (e% — ez)/(e.\ — e^). The function p(z) is periodic with respect to the two periods, 
2to and 2u/, as follows |29J30(1: 



p(z + 2Mu> + 2Nu') = p(z), 



K(k) 
V e i ~ e 3 



, ^ 



K'(k) 

V e l ~ e 3 : 



(39) 



for integer M and N. We may set (ei — 63) = 1, which yields 

,,2 _l r /C 

3^2 



3 7 2 + /C 

e l = ^^ > e 2 



6 7 2 
With these substitutions, equation ( |37| ) becomes: 



(37 2 ~/C) 
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(40) 
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P'(k) = 



X P (u) = 0, 

e -2/c 



r 



(41) 



This is now in the form of a Lame equation. We may solve it explicitly by introducing two dimensionless constants, 
a and b, by means of the transcendental relations: 



3p(a)+3p(b) = -p 2 , 

p'(a) = -p'(b), 



(42) 



where primes in this section will denote derivatives with respect to z = w/V e i — e 3 = u - Using the differential 
equation |29 30] 



p' 2 (z)=4p 3 (z)-g 2 p{z)-g 3 , 



g 2 = 2 (ex + e 2 + e 3 ) , g 3 = 4eie 2 e 3 , 



(43) 



the relations in equation ( [42| ) imply: 
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With these relations, solutions of equation (Ell) may be written: 



9 7 4 



<H(u + a/) ct(o/ + a) ct(u/ + 0) 



(44) 



(45) 



where a(z) and £(z) are the quasi-periodic Weierstrass functions, defined by the relations |2!j,|3(|: C( z ) — ~p( z ) an d 
a / (z)/a(z) = C,{z). The solution U p (u) is normalized as U p (0) = 1, and the linearly-independent solution is U p (—u). 
Using the initial conditions of equation (122), the fluctuations Xl( T ) can be written as a linear combination of U p (u) 
and U p (— u): 



Xi{u) 



1 



2^/2^(0) 



I i^L) U P (-u) + (l - i^gj- ) J7 P („) 



7^(0) 



7^(0) 



(46) 



where the primes here denote d/du — r y^ 1 d/dr. From the form of the number operator, equation (pij), solutions 
U p (±u) which grow in time will contribute to particle production. 

In general, solutions to second-order differential equations with periodic coefficients will obey Floquet's theorem pl| ; 
that is, for a periodic "potential" with period 2u (as is the case for equation (f4l|)), solutions behave as 



U p {u + 2uj) = U p (u)e lF{p \ 



(47) 



where the Floquet index F{p) is independent of time. The solutions U p (—u) have the Floquet index —F(p). The 
quasi-periodicity of the cr(z) functions pQ|, 



<j{z + 2oj) = — a{z) exp [2£(w) (z + id)} , 



(48) 



yields, for the solutions in equation 



F(p) = 2i [« (C(o) + C(6)) - (a + b) CM] 



(49) 



The solutions U p (u) will reveal exponential instabilities, then, whenever -F(p) has a non-zero imaginary component. 
Equations (g4|) and ( fhf ) reveal that it is these instabilities which are responsible for the highly efficient, resonant 
particle production during preheating. As is clear from the relations of equation (|42|), together with series expansions 

'(b) are real. From equation 



of ip(z), p'{z), and ((z) (see 1 30 

f 



such instabilities will occur only when 6oi/i p(6) and 



we see that for the physical wavenumbers < k < oo, there exists only one band in which these exponential 



instabilities may occur: 
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or, in terms of ip and the dimensionless wavenumber I: 

3 



Vo + 2/C^ + -K? + 2/C < £ 2 < y/Zipi + 6/C^ + 4/C 2 + 2/C. 



(50) 



(51) 



The existence of only one single resonance band for physical wavenumbers in this case matches the results found for 
Minkowski spacetime in |9j. This resonance band is in terms of comoving wavenumber k; because we have neglected 
interactions amongst the final-state bosons, the only way a mode Xi( T ) w ^ slide out of the resonance band in this 
system is when the backreaction AS(r) grows so large that it damps the oscillations of the zero-mode, halting the 
parametric resonance. We will demonstrate in section V that for preheating in an expanding open FRW universe, 
this range only includes subcurvature modes. 

In order to evaluate quantitatively values of the solutions U p (±u) in this resonance band, it is helpful to rewrite 
the solutions in terms of the quasi-periodic Jacobian theta functions, which possess very rapidly-converging series 
expansions. Using the relations p9| 



a(z) 



the solutions U p (u) may be written: 

U p (±v) = 



2^ 



i'&^v) exp 



C(")z 2 
2w 



ITT \ V + 



4^ 



i94(u + a)fl 4 (u + /?)0!(O) 

tiftv) M<») MP) 



exp [±ivF(p)] , 



where we have used luiv = u, Iujol = a, and 2w/3 = b. The relation |2£ 

w 2ojWi(w) 

further allows us to rewrite the Floquet index F(p) as 

Wa) , MP) 



F(p) 



M<») MP) 



(52) 



(53) 



(54) 



(55) 



With the constraints on a and b in equation ([42j) , the Floquet index will reach a maximum in the resonance band for 
a — > b. We may evaluate the Floquet index at this maximum by means of the series expansion for $i(v) |29[: 



Mv) 
M^ 



7TCtn(7Tu) + 47T 






.,2/1 



- (1 - q 2n ) 



sin(2mr?;) 



= 7rctn(7rt;) + 47rg 2 sin(27r«) + 0(q 4 



(56) 



where the elliptic nome q = exp[-iK'(k)/K(k)] satisfies < q < 0.0432. Writing F(p) ~ 2i(tf' 1 (/3 max )/tfi(/3 max )) near 
its maximum in the resonance band yields 



n< I ?^H ) = 4^ + 0^ 



Re 



Wmax) 



= -l6lT 3 q + 0(q 3 ), 



(57) 



where the maximum of the resonance occurs at /3 ma x- From the relations in equations (|42j) and ([14]), /3 max corresponds 
to the wavenumber 






(58) 



In other words, the maximum resonance will occur for quanta with wavenumbers near the top of the resonance band. 
This behavior also matches that found for Minkowski spacetime in ||. 

With these expressions, we may determine the growth of the backreaction due to created quanta, as well as the 
number of quanta produced during the preheating epoch. From the sign of F(p) at its maximum, it is clear that the 
Up(—u) modes will grow exponentially in time for modes within the resonance band, while the resonant U p (u) modes 
will decrease exponentially. Keeping only the growing modes, then, we may approximate 



\Xi(r)\' 



1 



8^(0) 



\U p (-u)\ 2 



1 



wftO) 



1 2 U' 2 (0) 



(59) 



and 



E(r) 



dit 2 .. a 

-^\Xi{r)\ 



dii 



16tt 2 w £ (0) 



\A p (u)Y 



<4(P) 



i 2 u p 2 (o) 



exp 



-Re(F(p)) 



(60) 



inside the resonance band. Here we have written the combination of oscillating ^-functions as the single function, 
A p (u). We will evaluate this integral by means of a saddle-point approximation. Using equation (J4l|), we may 
substitute 



dei^^dppJp*-^, 



(61) 



and using equation ( |57| ) we may set 

Re (F(p)) ~ 8nq - 16n 3 q (p - p max ) 2 + 0(q 3 ). (62) 

Furthermore, the integral will reach its maximum value when the oscillating term A p {u) = 1. This yields 



S(r) 



r'pjp 1 



■IK 



128tt 3 ^(0) 
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7 2 ^ 2 (0) 
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■. exp 
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(63) 



Following H , we may then write the backreaction in the form 

v 3/2 



E(r) 
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NilQyfF 



exp [B{K)^t] 



(64) 



In this form, we may solve for the time r when the backreaction 3AS(r) grows to be of the same magnitude as the 
tree-level terms. In our case, this occurs when 3AE(r) ~ 3?/> 2 (t). 

For the entire range — 2/C < "0? < oo, the average of the square of the zero- mode over a period of its oscillations 
gives, to a good approximation 



3V^C0^|^ = |(7 2 -/C) 



(65) 
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Setting the backreaction equal to this quantity yields 



7"cnd 



I / N(K.) ( 1 - 3/C 

■ In 



B{JC)l 



(66) 



2\JB(K) 



where we have defined K. = /C/(37 2 ). Once the backreaction grows to equal the magnitude of the tree-level terms, the 
parametric amplification of the fluctuation modes ends. This is the end of the preheating epoch, and hence we label 
the time at which this occurs r en d- 

In order to evaluate S(r), we must calculate U l2 (Q) at p max . Using the relations [p0| : 



c(. 1 + . 2 )=c(, 1 ) + c(, 2 ) + i- p,(zi) - p ' (Z2) 



2 p(zi)-p(z 2 ) 



then equations (|40| ) and (|45|) yield 

^; 2 (o)ip„ 



pV) = o . p( w ') = e 3, 
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[p{K 



e 3j 
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2 3 + 9/C ; 



3/2 



- 18JC 1 - K 



(3 - 3K\ Vz~- 



9/C 2 - 9/C 1 - K 



Similarly, the quantity 



W L x (°)=^ma X + 3^=7 2 

The backreaction may thus be written: 
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where we have defined 



G{t) 



j{£) = 



3 + 9K, 2 - 61C\/3 + 9JC 2 

3 + ^3 + 9K 2 - 3£ 
3 + V3 + QIC 2 ( 1 + iK. + 6IC 2 ) -itU-'dt- 2K, 



3/2 



(3 + 9/C 2 ) ' - 9/C (1 - A ■- 
When if = 0, this reduces to the more simple form: 

" 3/2 ^7J^t^± 



S(r)|K=o 



V"o 



: exp 



87i-y> q max 
^(l/v 7 ^) 



(67) 



(68) 



(69) 



(70) 



(71) 



2567T 3 ' y g max 

The numerical values K(l/y/2) = 1.854, q m3X = 0.0432, and A = 10~ 12 yield for the quantities in equation (p4|), 

256tt 3 



(72) 



AT(0) = 
B(0) = 



(In 



V3 + V3 V ^(VV2) 
87rg n 



K(l/V2) 



= 0.586, 



'T'end 



(0) 



■In 



AT(0) 



= 556.995, 



57.266 



B(0)Vo \2Av/B(0y 



V'o 



(73) 
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The quantities N(JC), B(JC), and r en d(^) will all approach these values in the limit tp — > oo. We will calculate the 
maximum quantitative deviation from the flat-space results in section V. 

Using the same saddle-point approximation, we may solve for the total number of particles produced during pre- 
heating. Within the resonance band, we will approximate 



dr 



—F(P) 

2uj yFJ 



Xt(r), 



so that near the center of the resonance band (see equation (J24|)), 



Nf(r) 



w*(t) 



\Xi(r)f 



.2„2„,21 



I6n 2 q 2 j 



.,2, ,,2 



(r) 



The total number of particles is 



Ctal(^) 



dee 2 

~2^ 2 



Nf{r). 



(74) 



(75) 



(76) 



From equation (p0|), using 3AS(r cn( j) — 3i/j 2 /2, and approximating w^(r en d) — <*^(0), this may be written 



N tatal( T end) 



(l - 3ic) 

4A 



3 + V3 + 9K 2 - 3/C 

i6ttV 
K 2 (k) (i + V3 + 9IC 2 - ilC 



1/2 



(77) 



Equation (177]) confirms that at the end of the preheating epoch, the number of particles produced is nonperturbatively 
large, iV tota i ex A -1 . Thus, for a massless inflaton decaying strictly into massless inflaton quanta, resonant preheating 
in an expanding FRW metric closely resembles the situation for Minkowski spacetime. 



V. COMPARISON OF OPEN AND FLAT CASES 



In this section, we will demonstrate that two key approximations made above are consistent for the case of preheating 
with a massless inflaton, and then compare the numerical results for the K = — 1 and K = cases. First consider 
the approximation made in section III that TL <C ijj during preheating. Taking 7i(0) = 0, we may write the average of 
this Hubble parameter over the duration of preheating. For K = 0, 



— 1 

it = ■^'HyTend) 



1 



•^Tcnd 



io- 2 ^o, 



(78) 



where the last step comes from using equation (73). As demonstrated below, for the K = — 1 case, r cn d remains 



greater than or equal to its numerical value in the K = case, so that 

H = -ctnh(r cnd ) < ip 



(79) 



for all ipo > — 2/C, as long as C 2 < O(10) after renormalization. Over the course of preheating, then, it remains a good 
approximation to neglect the Ti 2 ifj 2 term relative to the ip 4 term in the energy density of the zero-mode, equation 
(f27|); when this is done, the energy density remains proportional to a~ 4 when m = 0. 

The other main assumption to check concerns the expansion of the fluctuating field \ m subcurvature modes only, 
when K = — 1. For the decay process to be resonant, we require that the period of the Floquet solutions remain less 
than the Hubble time. In terms of the time-like variable u, this period is 2u — 2K(k); in terms of conformal time 77, 
then, this requirement becomes: 



2K(k) 



<H~ 



(80) 
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When K = — 1, a = V%H 1 , so this requirement becomes 

jC > 2V2K(k). (81) 

The modulus k is related to 7 as (see equation (J36[)) 



-2 t 2 C 2 - K 

k = V' (82) 

Equation ( |Sl| ) is therefore satisfied for 7C > 5.3. From equations ( ffl| ) and (|5C|), and using fc = C£, this means that 



3 / K 2 C 2 3 

fcmin = 2 V 7 4 C 4 + —3— + 2K > - (5.3) 2 - 2 = 40.1 , or fc min > 6.3. (83) 

Thus, over the entire range of allowable initial conditions for -0 O , preheating with K — — 1 will only populate subcur- 
vature modes, with fc m ; n > 1. 

With this restriction on the combination 7C, we may further calculate the quantitative deviation for K = — 1 in the 
growth of backreaction and total number of particles produced from the K = case. Numerical differences from the 
K = case arise both from the dependence of k (and hence of K{k) and q) on /C, as well as from the explicit factors 
of K, — K/(3^f 2 C 2 ) in such quantities as N, B, r cn( j, and N?*i ,. At the minimum allowable value of jC — 5.3, we have 
k = 0.72 (near the flat case of k = 0.707), or q = 0.0396 (near the flat case of q — 0.0432). With the constants at 
these values, 



From the definition of 7 2 , we also have the relation 



< 3^2 = 1-2 x 10- 2 . (84) 



^(l-m)- 1 ' 2 . (85) 

This means that the deviation from the K = case remains on the order of 10~ 2 . For example, comparing equations 
© and @, 



^maxiA^) 


= G(0) + 2.1 x 10" 2 


'J n\&,x_\i^' ) 


= J(0) + 1.0x 10" 2 


A r max(/C) 

N(0) 


= 0.93, 


-Bmax(^) 


= 0.89. 



B(0) (86) 

With these values, we see that the duration of preheating also remains nearly the same in the K = — 1 and K = 
cases: 

^§ = 1-14. (87) 

Tend(0) 

Finally, the total number of particles produced in the two cases obeys 

^t a o d t al (£) 



^tallO) 



= 0.98. (88) 



For 7C > (7C) min = 5.3, k — > l/-\/2 and K. ~ ■» 0, so all numerical quantities further approach their flat-space values. 

Thus, the requirement that the period of the Floquet solutions, 2oj, remain less than one Hubble time greatly 
restricts the numerical deviations between quantities in the K = — 1 and K = cases, when we consider chaotic 
inflation initial conditions. For all dynamically-consistent initial values of ip in the open inflation case, the total 
number of particles produced during preheating remains within 2 percent of the flat-space results. 
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VI. CONCLUSIONS 

For a massless inflaton decaying resonantly into massless infiaton quanta, preheating in expanding FRW spacctimcs 
closely resembles the Minkowski spacetime case. With chaotic inflation initial conditions, furthermore, quantitative 
differences in the spectra of produced particles remain small when comparing spatially-flat with open models of in- 
flation. In both cases, a single narrow resonance band for comoving wavenumber k exists. Preheating ends at T cn( j, 
when the backreaction AS(r) due to produced quanta damps the oscillations of the zero- mode. Because the Hartree 
factorization employed here neglects interactions amongst the final-state bosons, such as re-scatterings between each 
other and with the zero-mode, the quantitative values for the spectra derived here likely represent an overestimate of 
their true magnitudes, though with A ~ 10 -12 , such re-scatterings should remain subdominant effects. Furthermore, 
analytic study of the dynamics is useful for comparing preheating in the K — and K = — 1 cases. 

The analytic study also reveals clearly the limits for preheating with a massive inflaton, when the expansion of the 
universe is taken into account. From the form of the equations of motion in terms of the rescaled fields and conformal 
time, it becomes easy to understand why numerical studies of preheating in an expanding flat-FRW spacetime reveal 
a lack of resonant decays for massive inflatons and quanta, for many values of their couplings. |7|,18|| A non-zero 
mass for either the inflaton or the quanta would break conformal invariance. For chaotic inflation initial conditions 
as considered here, solutions for the zero-mode of a massive inflaton would no longer be simply-periodic or elliptic 
functions. Furthermore, if the mass of the zero-mode were zero but that of the decay products nonzero, then the 
equation of motion for the fluctuating field would include the non-periodic term, a 2 {rf)m 2 . Either of these situations 
would mean that the "potential" for the fluctuating field would no longer be periodic. Yet the existence of Floquet 
solutions, with one or more bands of exponential instabilities, depends upon a periodic potential in the equation of 
motion of the fluctuating field. Only in the limit gip 2 ^$> a 2 (rf)m 2 , where m is the mass of the produced quanta 
and g is the coupling strength between the inflaton and fluctuating field, will the potential approximate a periodic 
form. For the case of inflaton decay into inflaton bosons, as studied here, limits both on A and on m v from observed 
cosmic microwave background anisotropics make it impossible to satisfy this limit, and no preheating may occur when 
TO<p j^ for either K = or K = — 1, at least for chaotic inflation initial conditions. 

Finally, as demonstrated here, the K = and K = — 1 cases agree in the limit ip » |/C|, which is not surprising 
considering the equations of motion in equations ( J34| ) and (|4l]) . For this reason, it would be interesting to compare 
the preheating scenario for a symmetry-breaking potential with new inflation initial conditions, vb — ► 0, for both the 
K = and K = — 1 cases. The question of a non-thermal restoration of symmetry with such a potential and initial 
conditions has been raised for preheating in both Minkowski and flat-FRW spacetimes |15| (though see also P,|10|l). 
Because the preheating dynamics for such a model in an expanding open universe differ from those in both of these 
spatially-flat cases, such a scenario warrants attention. This is the subject of further research. 
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